The nature of the ordered phase of the confined self-assembled rigid rod model 
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We investigate the nature of the ordered phase and the orientational correlations between adjacent 
layers of the confined three-dimensional self-assembled rigid rod model, on the cubic lattice. We 
find that the ordered phase at finite temperatures becomes uniaxial in the thermodynamic limit, 
by contrast to the ground state (partial) order where the orientation of the uncorrelated layers is 
perpendicular to one of the three lattice directions. The increase of the orientational correlation 
between layers as the number of layers increases suggests that the unconfined model may also exhibit 
uniaxial ordering at finite temperatures. 



I. INTRODUCTION 



State of the art techniques for synthesizing colloids 
monodisperse in shape and size allow their collective be- 
havior to be investigated [l| . The new particles may be- 
come the 'molecules' of new materials if they can be tai- 
lored and assembled into useful structures Q. In fact, the 
possibility of particle decoration (through, e.g., glancing 
angle deposition, tcmplating, or lithography) produces 
colloids with pre-determined surface patterns (patches). 
Patches yield new features such as anisotropic interac- 
tions, control of the valence, and the formation of per- 
manent electrical dipoles, paving the way for the ratio- 
nal development of novel self-assembled materials (e.g., 
super-molecules) with highly tunable mechanical, opti- 
cal, and thermal properties (JQ. 

Self-assembly has been exploited theoretically for a 
primitive model of patchy colloids and state of the art 
simulation studies revealed how the number, type, and 
distribution of the patches determine the self-assembled 
structures. In systems with two bonding sites per par- 
ticle, only (polydisperse) linear chains form and there is 
no liquid- vapor phase transition (3|. If the linear chains 
are stiff they will undergo an ordering transition, at 
fixed concentration, as the temperature decreases. The 
minimal model of this transition considers the effects 
of the equilibrium polydispcrsity and the polymeriza- 
tion process of the rods. In this context, we proposed 
a model of self-assembled rigid rods (SARR), composed 
of monomers with two bonding sites that polymerize re- 
versibly into polydisperse chains Q and carried out ex- 
tensive Monte Carlo simulations to investigate the nature 
of the ordering transition on the square and triangular 
lattices d, @ . The polydisperse rods undergo a continu- 
ous ordering transition that was found to be in the two- 



dimensional (2D) Potts q=2 (Ising) and q=3 universal- 
ity classes, respectively, as in similar models where the 
rods are monodisperse 0- These findings refute previ- 
ous claims, based on Canonical Monte Carlo simulations, 
that equilibrium polydispersity and the statistical ensem- 
ble change the criticaiity of these models to random per- 
colation @-03 • 

The nature of the ordering transition of the three- 
dimensional (3D) SARR model on the simple cubic lat- 
tice is much more difficult to establish. The model con- 
sists of particles with two patches aligned along ±d, 
where a represents one of the three lattice directions 
(x,y,z). Particles on nearest-neighbor (NN) lattice sites 
Vi and Vi + a interact attractively with energy — e if their 
patches are aligned along a. Monte Carlo Simulations 
using efficient algorithms suggest that the ordered phase 
(below the transition temperature) exhibits a bias to- 
wards uniaxial behavior (i.e. the system exhibits a ten- 
dency to align different layers, by contrast to the ground 
state partial order). This tendency is observed only when 
the system size, defined by L with L 3 the number of sites 
considered in the simulation, is sufficiently large at tem- 
peratures that are not too low, T >> 0. Despite the 
use of efficient cluster algorithms we have not been able 
to establish the nature of the ordered phase as the sys- 
tem sizes required to observe uniaxial behavior increase 
rapidly as the temperature decreases, as discussed below. 

Here we consider the confined 3D SARR model as a 
first step towards elucidating the nature of the ordered 
transition on the cubic lattice. We investigate the nature 
of the ordered phase and the orientational correlations 
between adjacent layers of the confined model, on the cu- 
bic lattice, and find that the ordered phase at finite tem- 
peratures becomes uniaxial in the thermodynamic limit, 
by contrast to the ground state order where the orien- 
tation of the uncorrelated layers is perpendicular to one 
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of the three lattice directions. In addition, we find that 
the orientational correlation between layers increases as 
the number of layers increases from two to three suggest- 
ing that the unconfined model may also exhibit uniaxial 
order at finite temperatures. 

The paper is arranged as follows: In section II we de- 
scribe the ground state and the simulation methods used 
to analyze the 3D SARR model at finite temperatures 
while in section III we present the simulation results for 
the ordering transition and the order parameters. In sec- 
tion IV we introduce the confined SARR model. In sec- 
tion V we present the simulation results for the order 
parameters and the correlations between adjacent layers, 
for models with two and three layers. We conclude in 
section VI with a discussion of the results. 



II. THREE-DIMENSIONAL SARR MODEL 

A. Ground state 

In the full lattice limit every site is occupied by one 
particle aligned in one of the three lattice directions. In 
the ground state, the SARR model exhibits partial or- 
der: One lattice direction (say, z) is suppressed, with the 
uncoupled layers aligned in one of the remaining lattice 
directions (x or y). The ground state potential energy is 
then U = —Ne, where N is the number of lattice sites, 
with degeneracy: 



wgs = 3x (2 L - 1). 



(1) 



The entropy per site vanishes in the thermodynamic 
limit. 



B. Simulation procedures 

It has been shown that the full-lattice 2D SARR model 
on the square lattice can be mapped on the 2D Ising 
model 0| . This mapping allows the use of cluster algo- 
rithms developed for Potts models [ill, E3 to enhance the 
efficiency of the Monte Carlo simulations. In what follows 
we describe how the Swendsen-Wang algorithm may be 
adapted to the SARR model on lattices where the map- 
ping does not exist. We recall that the 3D SARR model 
on the cubic lattice or the 2D SARR model on the trian- 
gular lattice cannot be mapped on Ising or Potts models 
[f|. We can, however, develop cluster algorithms based 
on the layer structure of ground state of the 3D SARR 
model and the mapping of the 2D SARR model on the 
square lattice. 



1. Cluster sampling 

The cluster algorithm samples at each MC step a sub- 
set of all sites as described next. One of the three lattice 



directions is chosen at random, say z. Then the sites 
oriented along z are blocked, i.e. their state is frozen 
during the MC step. Sites with orientations x or y are 
active, and their states may change during the MC step. 
Given the NN character of the Hamiltonian the proce- 
dure may be (and it is) applied to all active sites in one 
MC step. For simplicity, however, we consider one layer, 
i.e. all sites, i, with Zi = zq. The procedure starts by 
checking the links between pairs of NN active sites (not 
to be confused with bonds of the original model). Two 
NN active sites with the same orientation are linked with 
probability B = 1 - exp(-/3e/2), where /3 = {k B T)- 1 . 
Links cannot be formed between pairs of NN active sites 
with different orientations. We define clusters of active 
sites based on the links generated in the previous step. 
The new configuration is obtained by choosing, indepen- 
dently, for each cluster a new in-plane orientation. The 
probability of the new cluster orientation, (x and y), is 
given by: 



Ak (a) cx exp 



-n k (a) 



(2) 



where (a) is the number of patches of the cluster k that 
point to blocked sites when the active sites are oriented 
along a. 

The procedure is validated using the plaquette 
formalism^ that maps the model with blocked sites ori- 
entated along £ to a Potts model [l~3|] in an external field. 
For a system with L layers (z = 1, 2, • ■ ■ , L) the intralayer 
potential energy can be written as: 



U 



<ij> 



5(6ti, a 3 )5{dii, fij); 



(3) 



where < ij > runs over the NN active sites in one layer, 
a.{ is the orientation of site i, and S(a, a') is one if \a-a'\ = 
1, and zero otherwise. This intralayer Hamiltonian can 
be mapped to a q = 2 Potts model[l3| with blocked sites 
on the square lattice. Using the plaquette formalism [|| 
we find: 



U = -K 8(&i,i 



<ij> 



Y^[Kv + Ki8{a u hk)]- (4) 

<ik] 



where the subindex < ik] runs over pairs of NN sites on 
the layer with i an active site and k a passive one. In Eq. 
(U]) K is the coupling constant, while Ko and K\ describe 
the interactions between active sites and the blocked ones 
(this may be viewed as the interaction of an external 
field with the active sites). K, Kq, and K\ are given in 
terms of the energy of the patchy model by: K = e/2, 
A'o = e/4, and K\ = —e/2. As expected the interac- 
tion energy of the active sites pointing to blocked ones is 
unfavorable. It is now straightforward to implement the 
cluster algorithm described above. The Potts coupling 
defines linking criterion between active sites according to 
the Swendsen-Wang rules fllj while the single-particle 
interactions are taken into account by considering the 
effect of an external field CDl . 



3 



2. Sublattice sampling 

In addition to the cluster moves we implemented sub- 
lattice single-particle moves. We consider systems with L 
even, and divide the sites into two sublattices: those with 
x i + Hi + z i odd (sublattice 1) and those with x± + yi + Zj 
even (sublattice 2). Notice that two NN sites belong to 
different sublattices. In one sublattice sampling move, 
we choose one of the sublattices at random and then 
update the state of each site by computing the inter- 
action with its NNs, Ui(a), for the three orientations: 
a = x, y, z. The new configuration is obtained by choos- 
ing, for each particle, a new orientation with probability: 
Pi(a) oc exp [—fiui(&)]. 

In order to check the cluster algorithm and its imple- 
mentation we have run pairs of simulations using cither 
cluster moves or sublattices moves only. The results were 
found to be the same within error bars. The cluster al- 
gorithm is much more efficient than the sublattice algo- 
rithm and the relative efficiency increases as the system 
size L increases. Nevertheless, as we will discuss later, 
its performance is far from optimal for very large sys- 
tems and temperatures slightly below the order-disorder 
transition. 

The Monte Carlo simulations of the 3D model are run 
in cycles. Wc choose at random, with equal probability, 
one of the five cycles to be run, namely (two) sublat- 
tice and (three) cluster samplings and then proceed as 
described above. 
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FIG. 1: (a) Excess heat capacities as a function of the tem- 
perature for different system sizes (three-dimensional SARR 
model), (b) Scaling of the peaks of the excess heat capacities 
with the system size. 



III. SIMULATION RESULTS FOR THE 3D 
SARR MODEL 



measures the uniaxial order: 



The order of the transition can be inferred from the 
scaling with the system size, of the peak of the excess 
heat capacity: c% x = (du/dT), where u = U/L 3 is the 
potential energy per site. At first order transitions the 
peak is expected to scale as: [l^] 



C~(L) = co + ^U 3 



4fc B T c 2 



(5) 



where T c is the transition temperature. In Figure [T] we 
plot the excess heat capacities in the transition region, 
and the scaling behavior of their peaks, c™ ax (L). The 
scaling of c" iax (L) with L 3 indicates that the transition 
is first order. A least-square fit yields the latent heat of 
the transition: AU/(Ne) = 0.024 ±0.001 . 

We consider two order parameters to characterize the 
transition. The first: O is based on the partial ground 
state order: 



= l-3min[N x ,N y ,N z ]/N, 



(6) 



where N a is the number of sites with orientation a, and 
N = L 3 is the total number of sites. The second: S 



S = 



Nl 



Ni 



N 2 



- 1 



(7) 



In the ground state O = 1 (one direction suppressed), 
while < S >= 1/4. Values of S greater than 1/4 be- 
low the order-disorder transition temperature signal the 
tendency for uniaxial order. In Figure [5] wc plot the two 
order parameters as a function of the temperature for sev- 
eral system sizes. Both sets of curves 0(T, L)and S(T, L) 
show, as L increases, an abrupt change at (or close to) 
the temperature where the heat capacity peaks. In small 
systems, S(T) varies monotonically with the tempera- 
ture. However, for the largest systems S(T) peaks at 
temperatures slightly below the transition temperature, 
and saturates at values larger than the ground state value 
at low temperatures. This finding suggests a surpris- 
ing tendency for uniaxial order, i.e. at finite subcritical 
temperatures and large system sizes the particles pre- 
fer to align in one direction rather than aligning in two 
directions as expected from the ground state analysis. 
Notice that as L increases the results for S(T) below 
T* = ksTc/e ~ 0.483 have large error bars. This is a 
signature of the lost of efficiency of the cluster algorithm 
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FIG. 2: (a) Order parameter O, and (b) Order parameter S 
for the three dimensional SARR model as a function of the 
temperature for several system sizes, L (see the legends). 



The confined SARR model exhibits an order-disorder 
transition as the models in 2D and 3D. Note that the 
limit of confinement (single layer) of the 3D model is not 
equivalent to the 2D SARR model as the patches can 
be aligned in three distinct directions. Given the nature 
of the ground state, where each layer is ordered in an 
arbitrary direction (x or y) we consider the single-layer 
order parameters Si, defined as: 



Si = -^[N x (i)-N v (i)} 



(8) 



where the index i refers to one layer. N x (i) and N y (i) 
are the number of sites on layer i with patches in the 
x and y directions, respectively. Due to the symmetry 
of the model, the average single-layer order parameters 
vanish, {Si) = 0. In the thermodynamic limit, at low 
temperatures, we expect an ordering transition described 
by: 



= 0; T > T c , 
> 0; T < T c . 



(9) 



We anticipate a discontinuous transition for systems with 
a large number of layers (as in 3D), and a continuous one 
for thin slabs (as in 2D). The continuous transition is 
expected to be in the 2D Ising class as in single layer 
systems. 



SIMULATION RESULTS FOR THE 
CONFINED SYSTEMS 



to sample S(T,L) slightly below T c as the system size 
grows. 

With the current algorithms and computational re- 
sources, however, we cannot establish the nature of the 
ordered phase at finite temperatures, in the thermody- 
namic limit. 



IV. CONFINED SARR MODEL 

In order to investigate the mechanism that may drive 
uniaxial order in the 3D SARR model, and to quantify 
it, we have considered a simpler model where the system 
consists of a number of layers (h) with L 2 sites. These 
layers are taken perpendicular to z direction. Periodic 
boundary conditions are only considered in directions x 
and y. The ground state of the confined SARR model 
is 2 /l -degenerate with all the particles in a given layer 
aligned along the x or the y direction. We note that 
geometrical confinement is also used to assist the self- 
assembly process of 3D systems and thus the study of 
confinement is also of some practical relevance [l4j . 

The confined systems are simulated using the algo- 
rithms described in section III Bl with minor adaptations: 
namely, the cluster moves are carried out only for layers 
perpendicular to z. 



In order to proceed we consider the scaling behavior of 
the heat capacity, or the related quantity u'p = (du/d/3), 
with u the potential energy per site. For systems in the 
2D Ising universality class, the scaling behavior at criti- 
cality is: 

^09 c )~lnL. (10) 

In addition we investigate the scaling behavior of the ra- 
tios g4i = (Sf)/(S 2 ) 2 , related to the Binder cummulants 
[l2| . For 2D Ising critical behavior g^{L,j5) for different 

L cross at the universal value g^ 1 ~ 1.168 [HI]. 

The results for one layer, h — 1 comply with the ex- 
pected 2D Ising critical behavior. The critical temper- 
ature, T c , is estimated from the Binder cummulant fol- 
lowing standard procedures Q. Considering system sizes 
in the range 12 < L < 96 we find T* = k B T c /e = 
0.5196 ±0.0001. This result is consistent with the behav- 
ior of the pscudocritical temperatures T C (L) defined by 
the peaks of the heat capacity as a function of L (results 
not shown). In addition, at the estimated T c the scal- 
ing of the average order parameter exhibits the expected 
Ising behavior: S 2 (L,T C ) oc L- 2 ^' (where f3' = 1/8 is the 
critical exponent for the magnetization). The results for 
<?4(L, T) and the scaling of S 2 (L, T) are shown in Fig. [3] 

For h = 2 and h = 3, the curves for different system 
sizes (12 < L < 48) cross at a value of 174, close to that 
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FIG. 3: (a) Scaling behavior of gi{T), and (b) Scaling behav- 
ior of the order parameter, S 2 for one-layer system systems 
(h = 1) and different system sizes. The dashed line in (a) 
marks the universal value g^f 1 for 2D Ising critical behavior. 
The crossings of the curves g 4 {T) and < S 2 (L,T) > L 1/4 for 
different system sizes confirm the expected 2D-Ising critical- 
ly. 



FIG. 4: (a) gu(T, L) for different system sizes, and (b) scaling 
of (du/d/3) for the confined model with h — 2. The system 
size dependence of the results suggest 2D Ising criticality (see 
the text). 



as: 



of the Ising universality class (See Fig0]). Somewhat sur- 
prisingly, for h = 3, the crossing of gu(L, T) of the inner 
layer occurs at a temperature slightly below that of the 
outer layers (plot not shown). This is likely to be a finite- 
size effect. A possible explanation is that, for these values 
of L, the different layers are almost independent; within 
this assumption the pseudo-critical temperature of each 
layer depends on the density of defects (number of sites 
with orientation z). The inner layer is expected to have a 
larger number of defects since these sites oriented along z 
can establish two bonds; by contrast, in the outer layers 
the sites oriented along z can form at most one bond. 
The larger density of defects reduces the stability of the 
ordered layer and thus its local pseudo- critical tempera- 
ture is lower. 

At finite temperatures, some particles will be aligned 
in the z direction. An interaction between adjacent layers 
results from bond formation between particles in different 
layers (z-bonds) and a correlation between Si and Si+i 
may appear. If present, these correlations may drive the 
bias to uniaxial behavior observed in the simulations of 
the 3D model. Let us define a global order parameter S 



1 " 



(id 



where h is the number of layers of the confined model. 
Again symmetry implies (S) = 0. The average value of 
S 2 may be written as: 



h h-l h 

E^ 2 > + 2 E E (Wj) 



(12) 



The correlation between two layers is defined as: 

(SiSj) 



[(Sf)(S 2 )] 



1/2- 



(13) 



The correlation depends both on L and T, Cij(L, T), and 
it is expected to vanish at low and high temperatures. 
Inspection of Eq. ([T2|) reveals that (S 2 ) behaves in the 
limit of low temperatures as: 



lim(,S 2 (L,T)} = ~. 

T->0 h 



(14) 
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FIG. 5: (a) Single layer order parameter for the confined 
model with h = 2. (b) Global order parameter for the con- 
fined model with h = 2. 
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FIG. 6: (a) Single layer (z=2), and (b) global order parame- 
ters for the confined model with h — 3 



We located the order-disorder transition of the model 
with h = 2 in the full lattice limit by considering the 
behavior of (du/d/3). In Figure 0] we plot the results for 
different system sizes. In the region around the maxi- 
mum we observe the scaling (du/d[3) ~ InL in line with 
2D Ising criticality. In Figure [5] we plot the results for 
the single-layer, (L,T), and the global, S 2 (L,T), order 
parameters for the same model. The single-layer order 
parameter exhibits the usual dependence on T and L. 
The global order parameter, however, exhibits a differ- 
ent behavior: the curves for different system sizes cross 
around T c , and then merge as the temperature decreases. 
The same qualitative behavior is observed for the con- 
fined model with h = 3 in Figure [5] While the scaling 
of (du/d/3) and gu suggests a continuous transition, the 
crossing of the curves < S 2 > (L,T) at criticality sug- 
gests a (weak) first-order transition, as the order parame- 
ter, \S\, in the thermodynamic limit, could exhibit a dis- 
continuity at the transition jumping from zero (T > T c ) 
to a finite value |5 C | > 0. 

In Figure [7^, we plot the correlation function c\2 be- 
tween the layers of the h = 2 system. As expected the 
correlation decreases and appears to vanish at low and 
high temperatures. The most relevant feature, however, 
is that the correlation between layers increases markedly 
with the system size, L. Figured reveals that the corre- 



lation increases as L 2 (in the range of sizes considered). 
These results suggest that for T > 0, in the thermody- 
namic limit L — > oo, the confined model becomes uniaxial 
(i.e., the layers will align along a unique direction). 

Now, we consider the effect of the number of layers on 
the correlation between adjacent layers. In Figure [8] we 
plot the layer-layer correlation Cij for h = 2 and h = 3, 
for systems with L = 32. Note that the correlation func- 
tions C12 and C23 are equal (except for statistical errors) 
due to the symmetry of the model. The main conclusion 
from the results of Fig [S] is that for a fixed value of L 
the correlation between adjacent layers increases with the 
number of layers. This suggests that the ordered phase 
of the three dimensional SARR model may become uni- 
axial, in the thermodynamic limit, at finite temperatures 
< T < T c . Note, however, that as the temperature 
decreases the system size required to observe uniaxial or- 
dering increases very rapidly. 

It is clear that the simulation algorithms used in this 
work loose efficiency as L increases at temperatures 
slightly below the critical temperature. In order to con- 
firm the trend to uniaxiality suggested by the results pre- 
sented so far, we return to the two layer system, and use 
an indirect method to compute the free energy difference 
A A = A xy — A xx , where the subscripts indicate configu- 
rations where the layers are oriented preferentially in the 
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FIG. 7: (a) Correlation function, C12 between the order pa- 
rameters of adjacent layers, h — 2. (b) Dependence of C12 on 
the lateral size of the systems with h = 2 at two temperatures. 
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FIG. 9: Free energy difference of the configurations of two 
layer systems where the layers are oriented in distinct and in 
the same directions. 



where AU(T) = U xy (T) - U XX {T). The potential ener- 
gies U a p (T) are computed using Monte Carlo simulation 
of relatively large systems, L = 64, L = 128 and L = 256, 
without cluster moves to avoid interconvcrsion between 
the two types of configurations. The results arc plotted 
in Figure IH1 A A increases with temperature and is pro- 
portional to L 2 . Thus for large systems and moderate 
temperatures, the confined SARR model is expected to 
exhibit uniaxial order. 




FIG. 8: Correlation functions between the order parameters 
of adjacent layers for h = 2 and h — 3, (with lateral system 
size L — 32). 

same (xx) and in different (xy) directions. AA is com- 
puted for large L using thermodynamic integration from 
low temperature (where the free energy of the two types 
of configurations is the same: AA(T ) = 0, as T — > 0). 
The free energy difference is then: 

1 J J T 1 



VI. DISCUSSION 

These surprising results may be interpreted as follows: 
Consider a two-layer system at low temperature with 
most particles aligned along the x or y directions. At 
T > 0, however, a number of particles will align along 
z. A z-bond lowers the energy by — e with respect to 
two independent z-sites, one in each layer, but isolated 
z-bonds do not contribute to the oricntational correla- 
tion between layers. Now, suppose that two z-bonds 
occur in NN positions (for instance one between sites 
r ai = 1) and r a 2 = (i,j,2), and the second between 
sites r bl = (i + 1, j, 1) and r b2 = (i + 2)) (See Fig- 
ure flU]) . At (low) subcritical temperatures this pair of 
NN z-bonds promotes the alignment of both layers in 
the direction defined by the pair (the x direction in this 
example as shown in Figure 1101 In practice, configura- 
tions with a different number of pairs of NN z— bonds 
along the x and y directions favor the alignment of the 
layers. A bond counting argument gives the probability 
of aligning one layer along the easy x direction over the 
probability of aligning it along y, when a single pair of 
NN z-bonds, along x, is present: 
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FIG. 10: Sketch of the effect of NN pairs of z-sites (or z- 
bonds) in the orientational correlation of the layers. Seg- 
ments represent particles with patches aligned in the plane of 
the layer, circles represent sites with patches aligned in the 
z direction. Crosses mark the in-plane bonds suppressed by 
the presence of 2— sites. Note that the number of such bonds 
depends on the alignment of the pair of z-sites with respect to 
the alignment within the layer. Therefore the configuration 
(a) has a lower energy than the configuration (b). 



The ratio of the probabilities of aligning the layers over 
the probability of not doing so is then: 



Pxx ~t~ Pyy 



cosh(/3e). 



(17) 



Pxy ~\~ Pyx 

An estimate of the density of NN z-bonds, a(T) is: 

a{(3) « e" 4 ^. (18) 



Note that only configurations where the number of NN 
z-bonds in the x and y directions are different contribute 
to the orientational correlation of the layers. Let us, 
however, consider the rough estimates given above. At 



T* = 0.40, pe = 2.5 and a(T) « 4.5 x 10~ 5 . For a sys- 
tem with L = 64 most configurations will not have NN 
2-bonds, and about one in five (0.186) will have one. An 
estimate of cyi is then: 



C12 a 0.186 x 



cosh(2.5) - 1 
cosh(2.5) + 1 



0.134. 



(19) 



which is close to the value obtained from the simulation 
C12 = 0.125 ± 0.007. This estimate supports the hypoth- 
esis that the uniaxial behavior results from the orienta- 
tional correlation between adjacent layers driven by the 
presence of NN z-bonds. 

The characterization of the ordering transition of the 
confined SARR model with h > 2 requires the develop- 
ment of more efficient cluster simulation algorithms and 
thus the behavior of the 3D SARR model cannot be inves- 
tigated at present. The problem is related, but not iden- 
tical, to the model for crystallization and vitrification of 
scmificxiblc living polymers investigated by Menon and 
co-workers in 2D and 3D [HEl- 
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